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We calculate the finite-temperature local spectral weight (LSW) of a Luttinger liquid with an 
"open" (hard wall) boundary. Close to the boundary the LSW exhibits characteristic oscillations 
indicative of spin-charge separation. The line shape of the LSW is also found to have a Fano-like 
asymmetry, a feature originating from the interplay between electron-electron interaction and 
scattering off the boundary. Our results can be used to predict how edges and impurities influence 
scanning tunneling microscopy (STM) of one-dimensional electron systems at low temperatures 
and voltage bias. Applications to STM on single-walled carbon nanotubes are discussed. 
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I. INTRODUCTION 

Metallic electrons confined to one dimension exhibit a plethora of intriguing effects, driven by interactions and 
the coupling to impurities and defects^. At low energies a clean system is described by the concept of a spinful 
Luttinger liquid (LL)^, with properties very different from those of a Fermi liquid: the quasiparticle pole vanishes 
identically and only collective modes remain, separately carrying spin and charge. The response of an LL to the 
addition of a local potential scatterer also differs dramatically from that of a Fermi liquid: The repulsive electron- 
electron interaction produces long-range density oscillations that get tangled up with the impurity potential in such 
a way as to suppress the single-electron spectral weight close to the impurity, as well as the conductance through 
i^^: In the zero-temperature limit and with a spin-rotational invariant interaction the impurity effectively cuts the 
system in two parts, with "open" boundaries (hard walls) replacing the impurity. The case of a magnetic impurity 
— which interacts dynamically with the conduction electrons — is similar: In the zero-temperature limit the physics 
is that of two LLs separated by open boundaries, with the finite-T response governed by a scaling operator that 
tunnels electrons through the boundaries^. The pictures that emerge in both cases are universal in the sense that all 
response functions depend only on the electron-electron interaction, with critical exponents which for a spin-rotational 
interaction is coded by the single LL charge parameter Kc- Details of the coupling of the electrons to the impurity, 
or the structure of the impurity potential, are irrelevant. 

The fact that an impurity in an LL drives the system to an open boundary fixed point^ has spurred considerable 
theoretical work on properties of LLs with an open boundary condition (OBC)i2424i2iiiiiSiiii4. Added interest comes 
from the fact that many measurements on one-dimensional electron structures — such as the single-wall carbon 
nanotubes (SWCNTs)'-^, or quantum wires, realized in gated semiconductor heterostructuresi^ or grown on metallic 
substratesii — are expected to be significantly influenced by electron scattering from the edges, where the confining 
potential to a first approximation can be treated as an OBC. 

Most work to date has focused on the local spectral weight (LSW) of an LL with an OBC, yielding predictions for 
single-electron tunneling and photoemission measurements close to an edgC"'^* or close to an impurity at sufficiently 
low temperaturesiS. Measuring the energy lo (with h=l) with respect to the Fermi level, the low-temperature LSW 
A{lo) close to an open boundary scales a.^/^^ 

A{UJ) ^ c^(^c-^-l)/2 (1) 

where Kc < 1 for a repulsive electron-electron interactionSS. This is to be compared with that of a clean system 
probed away from its edges, where A{lo) ^ 0;'^'^"'"^= )/4-i/22._ Experiments on SWCNTs seem to agree with the 
theoretical prediction that the tunneling rate of electrons should follow a characteristic power law with temperature^^, 
with a significant reduction of tunneling into the end of a tube as compared to tunneling into its interior ("bulk" 
regime Oscillation patterns that suggest spin-charge separation have also been seen in the tunneling conductance 
between two quantum wires produced by cleaved edge overgrowth^^, in qualitative agreement with theoretical results. 
In another line of research, photoemission spectroscopy measurements on quasi-one-dimensional organic conductors 
have been interpreted within a picture where the one-dimensional chains in the samples are cut by impurities into 
disconnected pieces, each modeled as an LL with OBCs. Again using results for the LSW, it has been argued^i^iiSiiLS^ 
that this approach gives better agreement with experiments than conventional theory where photoemission spectra 
are compared to predictions from ordinary "bulk" LL theory^^. However, this alternative interpretation remains 
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controversial and the issue has been difficult to settle, much due to the fact that photoemission measurements on 
these materials are subject to a variety of subtle effects. 

The most direct way to probe an LSW is via scanning tunneling microscopy (STM)iS. These experiments are deli- 
cate, as the STM tip must be positioned at a very small distance from the sample for electrons to tunnel^^ . While this 
is feasible for SWCNTs, the high-precision STM experiments that have been carried out have probed tubes deposited 
on metallic substrates. This leads to a suppression of the electron-electron interaction from screening charges, and 
early results were successfully interpreted within a free electron modetfS*^. In another effort STM measurements were 
performed on SWCNTs freely suspended over a trenchS^, thus bypassing the problem with screening charges. How- 
ever, the resolution achieved in this experiment was not sufficient to test for the expected LL scaling at small energies. 
In more recent experiments SWCNTs deposited on atomically clean Au(lll) surfaces were studied by high-resolution 
STM spectroscopy^, revealing that the electronic standing waves close to the end of a tube have en enhanced charge 
velocity which may imply spin-charge separation, and a fortiori LL behavior—. 

Turning to theory, the spectral properties of LLs with OBCs are by now fairly well understood, although some 
open problems remain. Maybe most pressing is the question about the very applicability of LL theory: What is the 
energy scale A below which the power law in Eq. becomes visible? Obviously, an answer to this question is 
essential for making sensible predictions for experiments. From numerical and other studies of the one-dimensional 
Hubbard model-^''- it is known that the decrease of the LSW — as predicted by LL theory — is often preceded by 
a sharp increase, and that this effect is particularly pronounced near an edgei^ti^ or close to an impuritjJ^S. The 
effect is expected to be generic for any one-dimensional metallic system where the amplitude for back-scattering is 
larger than for forward scattering. For some systems with a (weakly screened) long-range interaction, such as the 
carbon nanotubes, back-scattering gets suppressed above a threshold temperature, and one expects the asymptotic 
LL scaling in Eq. ^ to be visible at accessible energy scales, as is also suggested by experimental^. More work is 
needed, though, to obtain a reliable estimate of the crossover scale A, given data from the underlying microscopic 
physics. 

We shall not address this issue here, but rather revisit the problem of determining the full coordinate- and tem- 
perature dependence of the LSW of one-dimensional interacting electrons with an OBC, assuming that the energy 
scale is sufficiently low for LL theory to be applicable. Knowing the detailed structure of the LSW is important for 
making predictions of future high-precision STM measurements of LL systems, of which the SWCNTs are presently 
the prime candidates'^^. In earlier works the zero-temperature properties of the LSWi^, as well as the finite temper- 
ature properties of the uniform part of the LSW (neglecting Friedel oscillations)^, have been reported. Here we treat 
the full problem at a finite temperature and exhibit the LSW for different choices of interaction strength and band 
filling. We shall find that close to an open boundary the line shape of the LSW has a marked asymmetry as a function 
of energy with respect to the Fermi level, a property that arises from the phases that appear in the single-electron 
Green's function, and which has not been examined before. The form of the asymmetry in the neighborhood of the 
Fermi level resembles a Fano line shape, a feature expected universally whenever a resonant state (like that induced 
by a magnetic impurity in an electron system) interferes with a non-resonant one^. As we shall see, the origin of the 
asymmetric line shape in the present case is very different, and is formed by an interplay between electron-electron 
interaction and scattering off the open boundary. The asymmetry is fairly robust against thermal effects, suggesting 
that Fano-like line shapes produced by the reflection of interacting electrons off boundaries can be observed at tem- 
peratures higher than those originating from their interference with a resonating level. Having access to the full LSW 
we will also be able to give a systematic description of how charge- and spin separation shows up as an oscillation 
pattern when being close to an edge (or, an impurity, at low temperatures). This information, which we extract for 
different temperatures, can be directly translated into a prediction of the measured differential conductance when 
probing an LL system by STM. Also, given the full LSW we derive its crossover from boundary to thermal scaling 
near the Fermi level. The thermal effects soften the power law singularities of the LSW, since the non-chiral terms 
in the electron Greens's function produce a leading scaling term that is linear in energy at higher temperatures. 
This softening should not be confused with the averaging effects that always occur when the experimental tunneling 
currents are calculated by integrating over the Fermi-Dirac distribution. 

Our paper is organized as follows: In Sec. II we review some basics about temperature-dependent local spectral 
weights and STM currents. In Sec. HI we derive an exact representation of the local spectral weight for a Luttinger 
liquid with an open boundary, paying due attention to the phase dependence that has not been examined in earlier 
studies. In this section we also show how to adapt the theory for applications to scanning tunneling microscopy of 
SWCNTs. Sec. IV contains our results, and in Sec. V we summarize the most important points. A reader mostly 
interested in the physics of the problem is advised to go directly to Sec. IV. Unless otherwise stated we use units 
where h = ks — c — 1. 
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II. PRELIMINARIES 



In order to calculate tunneling currents e.g. from an STM tip we will consider the transition rate of adding electrons 
to an LL system at a position x and with energy w, 

r+(L^,a:;/3) = 2VZ-i^exp(-/3^„)|(n|^'t(x)|m)|2^(a;-S„+i;„). (2) 

This expression follows from Fermi's golden rule assuming a tunneling Hamiltonian of the form — (7\E'J.-0o-,tip + h.c, and 
treating the tip as a reservoir with unit probability that an electron is available for tunneling. Here Z is the partition 
function of the A'^-particle system, \E'J.(a;) creates an electron in the sample with spin a at a:, and "00- tip removes an 
electron of the same spin from the tip. Equation (01 represents the probability that the A^-particle states | m) of energy 
Em are connected to the {N + l)-particle states | n) of energy = + w by the addition of an extra electron of 
energy to and coordinate x. The transition rate T~ of removing an electron is given by Eq. |2Jl by simply replacing 
the index m by n in the Boltzmann weight, assuming that a "hole" is available in the tip with unit probability. 

In order to calculate the transition rates it is useful to define the single-electron local spectral weight (LSW) 
A{uj,x]f3) which is directly related to the transition rates by 

A{uj,x- /3) = Z-\l + exp(-/?c^)) ^exp(-/3S,„) (a,)|™)|2 ^(^ _ + e„,) 

= {l + exp{-Puj))T+{uj,x;l3)/2TTg^ (3) 
= {l + exp{(3iu))r-{uj,x;f3)/2TTg\ 

It is well-known that the LSW defined in this way can be extracted from the spectral representation of the single- 
electron retarded Green's function 

G"it,x,P) = -ie{t){{^,{t,x),¥^{0,x)})f3, (4) 

by using tha^ 

A{uj, x; (3) = Im / G^(i, x; (3)e"^^dt. (5) 

At zero temperature this quantity is known to be the single-electron local density of states N{uj, x) in agreement with 
the definition in Eq. (Q. 

We shall extract the LSW in the standard way by first calculating the single-electron retarded Green's function. 
The calculation of A(uj, x; P) for the present problem requires some care in order to analyze the analytic structure 
and phase dependence in detail. In fact, the result which we derive in the next section, using bosonization, reveals 
a surprising asymmetric energy dependence of the LSW close to the Fermi level for a semi-infinite LL with an open 
boundary condition (OBC). 

Before taking on this task, let us recall how scanning tunneling microscopy (STM) is used to experimentally probe 
the LSW close to edges and impurities. In the simplest approach, when the STM tip is assumed to couple only to 
the conduction electrons (thus neglecting tunneling into localized impurity levels) the tunneling current is given by 
the integrated difference between the transition rates r+ and F^, weighted by the corresponding probabilities that 
an electron [hole] is available in the tip for tunneling to [from] the sample. With an applied voltage V one thus has: 



IiV,x;l3) = e NsTAiiu^ - eV) [/{lo - eV)T+iu;,x;p) ^ {1 - f{uj - eV))T-iu;,x;P)] da; 

'J —OO 

poo 

27Teg^po [f{LU~eV)- f{uj)]A{uj,x;f])dij, (6) 



where /(w) is the Fermi-Dirac distribution and we have approximated the density of states NsTAiii^) in the tip by a 
constant po in the last step. It is clear that we recover the conventional formula for tunneling at zero temperature^ 

I{V, x) = 2eng^ / iV(w, x)Nstm{uj - eV)duj, (7) 
Jo 

where N{oj,x) is the local single-electron density of states for a conduction electron in the sample, and iVsTM(w) is 
the density of states of the STM tip measured relative to the Fermi energy. By differentiating, the local differential 
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tunneling conductance can then be directly related to the local density of states in Eq. Q 

^2e'Tig^p^N{V,x). (8) 



dI(V,x) „ 2_„2, 



dV 

This expression remains valid at a finite temperature T, provided that the thermal length Ay ^ Vg/T is larger than 
any other characteristic length L of the experimental setup (such as the distance between the STM tip and the edge 
of the sample). The speed Vg that determines At is that of the spin collective modes (which in a one-dimensional 
interacting electron system are slower than the collective charge modes). When L > At a temperature-dependent 
description becomes necessary, and the expression for I{V,x) has to be modified according to Eq. The local 
differential conductance at finite temperature can therefore be written as 

dm /3) _ 2^^2^^^ r eV)A{uo, x; P)d^. (9) 



dV '^"J-oodV 

It follows that the line shape properties of the local tunneling conductance are directly determined by the LSW. 

With these preliminaries we now turn to the calculation of the finite-temperature LSW for an LL with an open 
boundary. 



III. DERIVING THE LOCAL SPECTRAL WEIGHT 



We consider an interacting electron liquid on a semi-infinite line, x > 0, subject to an OBC at the end x — 0. 
Following standard Luttinger-liquid approach^, we linearize the spectrum and decompose the electron field into 
left- (L) and right- (i?) moving chiral fermions at the two Fermi points zLkp, 

^,(x) = e-^^-^^j^Jx) + e^'^-'^^jj.Jx). (10) 

The zero-temperature single-electron Green's function at a point x can then be expressed in terms of the propagators 
of the time-evolved chiral fermions 

Git > 0,x) = {^,{t,x)^UO,x)) = (V'L.(i,x)vL(0,x)) + (V'^<,(t,a:)^lj,(0,x)) 

+ e*''^^(^i^.(^,^)V'L(0,x)) +e-'2'^---(V'^,(i,x)^tj^(0,a;)). (11) 

We see that there are two types of contributions to G{t > 0,x): oscillatory and non-oscillatory. While the latter 
are always present, the former are nonzero only if the left- and right moving fermions get entangled at a boundary. 
Imposing an open (Dirichlet) boundary condition at the "phantom site" which is situated one lattice spacing a from 
the end of the LL at a; = —a, 

M,4-a) = e^'=->^,(-a) + e^^'^-^^P^J-a) = 0, (12) 
and assuming that the chiral fermions are slowly varying on the scale of a, it follows that 

^rAO) = e^'^^M (13) 

where 

7 = TT + 2fc_Fa = 7r(l + ne), (14) 

with Tie the filling factor {ue — l for a half-filled band). Although not essential here, the "softening" of the boundary 
— implied by imposing the Dirichlet condition at a; = —a — is sometimes useful for modeling the dependence of the 
scattering phase shift 7 on the shape of the edge- or impurity-potential. The value of 7 may therefore depend on 
the details of the boundary geometry, but it is important to notice that it is in general not a multiple of n even at 
half- filling. Using Eq. (|13|l to analytically continue to negative coordinates^, the right-movers may be represented by 
left-movers as 

^jiJx)=e''^i^Lj-x), x>0. (15) 

We can then express the Green's function in Eq. 111|) in terms of left-moving fermions only, now taking values on the 
full line —00 < X < 00 

G{t > 0,x) = (VL,(i,x)^L(0:^)> + (V'L,(i, -2^)^-1,(0, -x)) 

+ e*(2'=.-+7)(^^^(i^ -x)^LiO, x)) + 6-^(2'=-+^' (^^,(i, a;)VL(0, -x)). (16) 
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Introducing 

GLL{t,x,x') = {^LAt^x)i:l^{0,x')) = (V'L(i,a:)Vi,(0,a.')>, (17) 

with the second equahty foUowing from the charge conjugation symmetry of the hnearized theory, Eq. ]l<j\i may be 
written as 

Git >0,x)^ GLLit, X, x) + GLL{t, -X, -x) + e'(2feF^+7)G'^^(t^ ^) + e-'<^^''^^+'>^GLLit, x, -x). (18) 
With the definition in Eq. Q the retarded Green's function can finaUy be cast on the compact form 

G"{t, x) = -iQ{t){AReGLL{t, x, x) + 2e'^^'"=-''+"''> ReGLL{t, -x, x) + 2e-''^^'""'+^'> ReGLLit, x, ~x)), (19) 

using Eqs. (|17(l and H18|l. To obtain the LSW in Eq. © we thus need to calculate the chiral Green's function in l|17(l . 
identify its real part, and then Fourier transform the resulting expression for G^{t,x) from H19|l . The first part can 
be done analytically by using bosonization, and we turn to this task in the next section. 

A. Chiral Green's function from bosonization 

Using standard bosonization^^ we write the left- and right-moving fermion fields as coherent superpositions of free 
bosonic charge and spin fields, (l>rc — {4'ri + 4'ri)/V2 and (/)rs = (0rt ~ 4'ri)/V2, with r = L,R 

V'Lcr(^'^) ^ -^1^^ exp (-iV2n{cosh9 4'Lc{x,t) + sinh6 (l)Rc{x,t) + a(j)Lsix,t))] 

= ^^^exp uV27r(cosh6'(/i_R,c(a;,t) +sinhe'(/)Lc(a;,t) -l-cr0fls(a;,t))) . (20) 
\I2tioi ^ ' 

Here a is a small-distance cutoff of the order of the lattice spacing of the underlying microscopic model, and r\ra are 
Klein factors obeying a diagonal Clifford algebra that ensure that fermion fields of different chirality r and/or spin 
cr anticommute. The parameter Q is related to the LL charge parameter Kf^ by = e^^, and is parameterized by 
the amplitudes for the low-energy scattering processes. For a system with long-range interaction these amplitudes 
become momentum dependent, but since we shall only be interested in the asymptotic low-energy behavior of the 
Green's functions we can restrict them to zero momentum and treat as a constant (taking a value < 1 for repulsive 
interaction). Note that this shortcut assumes a finite-range interaction, whereas an unscreened Coulomb interaction 
which diverges at vanishing momentum leads to very different physical. Away from a half filling umklapp scattering 
vanishes, and standard RG arguments show that backscattering processes become irrelevant (again assuming a re- 
pulsive electron-electron interaction). As is well-known, for this case the remaining dispersive and forward scattering 
vertices can be written as quadratic forms in bosonic operators, leading to two free boson theories, one for charge, 
and one for spin 

^ = E Y ((^-'^i.f + iP^Rif) ■ (21) 

J=c,s 

This defines the LL Hamiltonian density, here expressed in the chiral fields, with Wc(s) the speed of the charge (spin) 
bosonic modes. 

The logic of the construction just sketched is strictly valid only for a translational invariant system where all 
interaction processes can be classified into dispersive, forward, backward, or umklapp scattering^. For a system 
with an open boundary, translational invariance is broken and a two-particle interaction leads to additional scattering 
processes. As shown by Meden et ali^, however, the theory in Eq. H21|) still captures the universal low-energy physics. 
Perturbative arguments suggest that the energy range where it applies increases with the range of the interaction of 
the original microscopic theory. 

To make progress we analytically continue the charge 4>lc and spin (j)Ls boson fields in (|20|) to a; < such that the 
boundary condition in Eq. (|13|l is satisfied 

(t)Lc{t,-x) = -<j)Rc{t,x) + 



V2ii?c 

(l)Ls{t,-x) = -~<j)R,{t,x), x>0. (22) 
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Using ^ we can write 



X exp i\/27r(cosh6' (/)Lc(i, x) — sinh ^? (/)Lc(i, —a;) + a(f>Ls{t, x))^ 



(23) 



Given l|23|l the chiral Green's functions in p9(l are now easily calculated, using that (j)Lc/s are chiral bosons governed 
by the free theory in ll^T|> ^ 
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(24) 
(25) 



with GLL{t, —x, —x) and GLL{t, —x, x) obtained from l|24|) and l|25|) respectively by taking x — > —a;. The exponents 
are given by 



ki ^ iK, + l/Kc + 2)/8 
k2 - (Xc + - 2) /8 

k3 = {l/Kc-Kc)/8 



(26) 



To obtain the finite-temperature Green's function we use conformal field theory techniques^ and map the complex 
planes {zj = VjT + ix} on which the zero-temperature chiral theory is defined (with T = it the Euclidean time, and 
J = c, s) onto two infinite cylinders Tj = {w = Vjt' -\- ix'} of circumference (3 = 1/T 



Employing the transformation rule 



arctan 



,3«j 



1 = 1 



dz 



(27) 



(28) 



we obtain for the finite-temperature versions of the Green's functions 
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(30) 



We have here dropped the primes on the transformed coordinates and reinserted the real time variable. 



B. An exact representation of the LSW 

Having obtained the chiral Green's functions we next calculate their real parts, which, according to Eq. p9|) . define 
the LSW. It is essential here to consistently identify the phases of the Green's functions that appear in (|19|l . This is 
most easily done by exponentiating the expressions in H29|l and (|30|) and choosing the negative real axis as branch cut 
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FIG. 1: [color online] Branch points and branch cuts for GLL(t,x,x; 0). The phases are different in different regions defined 
by branch points. 



of the logarithm'*^. Treating i as a complex variable, this amounts to choosing the branch cuts for GLLit,x,x; (]) as 
shown on Fig.^ Functions GLL{t,x,x; (3) and GLL{t, x, —x; P) have different phases in different regions defined by 
branch points, since the phases differ by 27r on opposite sides of the cuts. Taking the cutoff a — > we get 
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with /ci, A:2, and defined in Eq. (|26|) . and where 



TTk2 

2 
7rfc2 
2 



0<t < 



at) = 



2 2 
+ ^ + ^ + TTfca 



Trfc. t < 



2x 



2x 



<t<- 

<t<0 



(33) 



Tvki 7rfc2 

2 2 



We see that the phases that appear in the real parts of the Green's functions take different values in different domains, 
consistent with the fact that the end points of the domains are branch points (see Fig.^). Given the results in Eq. 
(I32|l we can finally write down an exact representation of the LSW for an LL with an open boundary. Combining 
and ifT^ we find that 



4 f°° 

A{uj,x,l3) = — dtcosujtReGLL{t,x,x; j3) 

2 f°° 

+ — dtcos{2kFX + — ujt)ReGLLit,x,—x; p). 

J-oo 



(34) 



This expression reveals that A{(jj, x, f3) has a nontrivial dependence on the interaction strength (via ki,k2, k^ and the 
velocities Vc and Vg that parameterize the chiral Green's functions), in addition to an oscillation in the second term 
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that is shifted by a phase controlled by the band filling via kp and 7 in Eq. (|14|l and the distance x to the boundary. 
This leads to an interesting asymmetric dependence of the LSW on energy as will be discussed in Sec. IV. 

One of the most promising candidates for LL behavior are single- walled carbon nanotubes (SWCNT), which however 
require some modification to the Green's functions in Eqs. H29l) - (|32|l due to the presence of two electronic channels. 
Depending on the substrate the interaction constant Kc may be very small ~ 0.25 and strongly k-dependcnt for 
isolating substrates or closer to unity « 0.6 — 1 for metallic substrates. Especially in the latter case backscattcring 
processes may also play an important role at very low temperatures. Following '42"43'| the low energy expression 
for the electron field is written in terms of Bloch functions on the two graphite sublattices in order to arrive at two 
spinful channels of ID Fermion operators ipLaa and tpRaa where a =1, | is the spin and a = ± labels the two distinct 
channels. In order to bosonize the problem it is then possible to introduce bosonic fields of total and relative {5 = ±) 
charge and spin chiral bosons (fi^jg (r = L, R, j = c, s) and write the bosonization formula 

IpLaa = r-— (^^P{~'iV'^[cOSh9(pLc+{x) + SlB-h 9 (j) iic+ (x) + a(t>Lc.-{x) + (J(t)Ls+{x) + aiJ(t)Ls-{x)]), (35) 

v27ra 

where ipRaa is obtained by taking the complex conjugate of the right-hand-side of Eq. H35|l and switching L ^ R. 

An open end or an impurity will in general mix the two channels so the effective "open" boundary condition will 
in most cases not be a simple rcfiection. The analytic continuation analogous to Eq. I|22|) may therefore turn out 
to be more complicated for the four bosons. In the simplest symmetric cases we see that the boundary Green's 
functions have again a structure as given in Eqs. (|29I32I) . however with the exponent 1/2 replaced by 3/4 and 
ki = {K,+ + l/K,+ + 2) /16, fc2 = {Kc+ + l/K,+ - 2) /16 and fcg = (l/ifc+ - K,+) /16. Accordingly C(t) and C'(t) 
are given by Eq. (|33|) with each occurrence of the term 7r/4 in H33|) replaced by 37r/8. The oscillating Friedel-like 
terms in the LSW in the second term of Eq. (|34|l turn out to have the opposite sign on the two sublattices^. 

IV. RESULTS 

In order to extract the physics from the LSWs derived in the previous section we have numerically carried out the 
Fourier transforms in Eq. (|34|l . The results reveal a rich structure in the LSW of a Luttinger liquid when an edge 
or an impurity is present. We here focus on three aspects of the LSW of particular interest for STM experiments: 
Its asymmetric line shape as a function of energy, oscillation patterns revealing spin-charge separation, and thermal 
effects. To keep the discussion general, we do not specify the dependence of the velocities Vc.s on the Luttinger liquid 
parameter Kc, which is sensitive to the choice of microscopic model. Since our results do not depend on the particular 
relation between Vc.s and Kc, we choose them freely so as to be able to observe all possible scenarios for the behavior 
of the LSW. The question whether a particular scenario will be realized in a particular model is of course determined 
by the relation between velocities and Luttinger liquid parameter. 

A. Properties of the LSW 

One of the most striking feature of the LSW of an LL is the well-known power law suppression at low energies 
proportional to lo'^^''^^'' For the special case of zero temperature and commensurate oscillations (i.e. setting 

2kpna -I- 7 = ttto (m S Z)) the LSW of Eq. H34II near a boundary at x = na has been analyzed befor o^i^^ . In that 
case the first term of Eq. (|34|) takes on a scaling form with the variable ujx that shows a crossover from bulk scaling 
to boundary scaling in Eq. |^ with slow oscillations proportional to siii{2ujx/vc If the phase is 

neglected or set to 2kpna -I- 7 = tott the second Friedel-like term in Eq. 134|) also takes on a scaling form, but shows 
both spin and charge modulations that decay with aj'^^^^x^'''^^^^^ . We will now focus on the energy dependence at 
both small and large temperatures together with the effect of the phase 7. 

1. Asymmetric line-shape 

At small energies, the most striking feature of the LSW is the asymmetry of its line shape, as seen in Fig. |21 From 
Eq. (|34|) it is apparent that this property is due to the phase 2kpx + 7 appearing in the second term of A(uj, x; (3). 
Here the distance from the boundary x = na has to be measured in integer units of the lattice spacing a, since the 
original Fermion operators typically correspond to orbitals in the crystal lattice. The phase causes a shift of the 
periodic structure of A{u!, x; (3) with respect to the Fermi level, and accordingly determines how the spectral weight 
suppression close to the Fermi level affects the line shape. The shift depends upon the filling factor Ue since, for fixed 
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X, Ue determines the values of fc^? and 7 in Eq. p4|l . The LSW is asymmetric in general, but for particular fillings 
for which 2kpna + 7 = 7rm (m g Z) it becomes a symmetric function of uj. By inspection of Eqs. H32|) and H34|) one 
finds that for a given interaction strength the asymmetry tends to zero very close to the boundary {2lox/vs <^ 1) 
and very far from the boundary {2lox/vc 3> 1, "bulk" regime) and reaches a maximum in the intermediate region. 
It is important to realize that the shift of the periodic structure with respect to the Fermi level is present also for 
the non-interacting case when the LSW takes the simple form A{uj,x](3) ~ cos{2kpx + 7 — 2ujx/vf)- We conclude 
that the shift is a pure boundary effect, and is due to the interference of the incoming and reflected electrons at the 
boundary. In contrast, the dip of the spectral weight at the Fermi level is an interaction effect. 

In Fig. 131 we show the energy and coordinate dependence of the LSW for even and odd sublattices. The Friedel 
oscillations on the scale of the lattice spacing a are easily visible as a flip of the asymmetry when going from one 
graph to the other. The spin and charge modulations of the amplitude are also present over longer wavelengths in real 
space, but are not clearly visible due to the relatively narrow coordinate range in Fig. |3| Note that the asymmetry 
with energy varies with the distance to the boundary since the phase 2kpx + 7 also varies with distance. 

It is interesting to note the similarity of the typical asymmetric line shapes in Figs. [3 and |31 with that of a Fano 
resonance-^^ when uj is close to the Fermi level. A Fano resonance is known to develop in the LSW for non-interacting 
electrons when coupled e.g. to a magnetic impurity, the effect being produced by the interference between resonating 
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FIG. 4: [color online] A{lj,x;P) for A'c = 0.7, v^/vf ~ 1.43, Uc/'fJs « 3, T/kpVF = 2.6 x 10"*, = 0.97. 



and nonresonating electron paths through the impurity-^ . As we have seen, the asymmetry in the present case instead 
comes from the combined effect of electron interactions (causing a dip in the LSW at the Fermi level) and the reflection 
of electrons off the boundary (causing a phase shifted oscillation in the LSW). As one expects the Fano line shape to 
survive for interacting electrons coupled to a magnetic impurity it is indeed satisfying to see this feature reproduced 
by the open boundary for which the impurity gets traded at low temperatures. 




-0.2 -0.1 0.1 0.2 -0.2 -0.1 0.1 0.2 



FIG. 5: [color online] A(lj,x;/3) for T/kpVF = 2.6 x 10"®, rie = 0.97, with Kc = 0.7, «c/wf ~ 1.43, Vc/vs « 3 (a) and = 0.9, 
Vc/vf ~ 1.11, Vc/vs ~ 1.26 (b). The color coding is the same as in Fig. 0| 



2. Spin-charge separation 

The proximity to an open boundary reveals a key property of interacting electrons in one dimension — spin-charge 
separation — i.e. the fact that the collective spin and charge excitations (induced e.g. by inserting an extra electron 
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into the system), propagate with different speeds, and hence "separate"—. The effect shows up in the LSW as a 
characteristic peak structure at intermediate distances from the boundary. Very far from the boundary {ujXg 3> 1) 
j4(r, w;/3) is a monotone function scahng as near the Fermi level, with bulk exponent a = {Kc + K^^)/A — 1/2 
at low temperatures-^. Extremely close to the boundary {ujXs <C 1) A{r,Lu;f3) has a similar structure, but with an 
enhanced suppression near the Fermi level, A{r,Lu;f3) ^ uj'^^ , with boundary exponent as — {K^^ — l)/2 at low 
temperatures^i^ (see Fig. 0] with color coding in arbitrary units). As one moves away from the immediate vicinity 
of the boundary an oscillation pattern emerges, which becomes most pronounced when ojXs ''^ C(l)- This oscillatory 
feature (see Fig. ISJ which is a superposition of spin and charge waves is due to the two types of branch points in 
Eqs. and which define the propagating collective spin and charge modes. The two panels correspond to two 
different choices of K^, with two different values of the ratio Vc/vs, leading to the two different peak structures. By 
close inspection of the graphs one can easily read off the corresponding velocity ratios. To see how, let us make a 
"cut" of the two panels in Fig. [S] at a distance x = 100a from the boundary, yielding the panels in Fig. El 





CO / kpVp 



CO / kpVp 



FIG. 6: A{Lo,x-f3) for T/kpVF = 2.6 x 10"^ = 0.97, with x = 50a, = 0.7, Vc/vf ~ 1.43, Vc/vs = Xc/Xs ~ 3 (a) and 
X = 100a, Kc = 0.9, Vc/vF « 1.11, v^/vs = Xc/Xs ~ 1.26 (b). 



When Kc — 0.7 (panel (a)) the spin and charge velocities differ significantly. For this case the short wavelength 
(As = TTVs/x) spin oscillations are modulated by long wavelength charge oscillations (Ac — ttVc/x). We note that 
there are three spin oscillations per one charge oscillation which is in agreement with the fact that in this case 
Xc/Xs = Vc/vg ~ 3 (see Fig. El (a)). When Kc gets closer to unity (non-interacting limit) the spin and charge velocities 
approach each other and the spin-charge separation manifests itself as a beating pattern, provided that Kc is not 
identical to unity (see Fig. El (b) where Kc = 0.9). The short wavelength (A = 27tVsVc/x{vs +Vc)) oscillations are now 
amplitude modulated by long wavelength oscillations (A' = 2'kVsVc/ x{vc — Vs)). We see that there are about six short 
wavelength oscillations per "bubble" of the long wavelength amplitude modulations, which is in agreement with the 
fact A'/A « 12 in this case, corresponding to Xc/Xg — Vc/vg ~ 1.26. 

The spin-charge oscillations in the LSW are also present in real space as discussed before^^. The oscillations as 
a function of energy might be easier to observe, however, since here there are no superimposed Friedel oscillations. 
We also remark that the existence of the beating pattern for values of Kc close to unity was proposed in Ref. ^3 
as a diagnostic tool for spin-charge separation in possible tunneling experiments of an LL where a scanning probe 
microscopy tip would be used as an impurity. 

One can map the dispersion of the spin and charge waves by taking a Fourier transform of the LSW (see Fig. [71 
(a)). This Fourier transform should not be confused with the momentum or angle resolved spectral weight, which is 
measured in photoemission experiments, although it shows similar features of spin charge separation. The dominant 
weights in the transform correspond to the uj(k) dependence of the excitations. The dispersion lines at fc = come 
from the non-oscillatory part of the LSW and represent the charge excitations, since the non-oscillatory part contains 
only charge oscillations. This feature at fc = is also a clear indication of interaction effects and disappears as Kc 1. 
The dispersion lines at fc 7^ come from the Friedel terms, and contain spin and charge branches shifted from k = 
by zL2kF- The mirror symmetry about fc = reflects the standing wave nature of the oscillations. In Fig. [71 (b) the 
dispersion relations uj{k) — ztkvc/2 at fc « and w(fc) — (±fc — 2kp)Vci s/'2 at fc « ±.2kp are plotted on top of the 
Fourier transform and agree with the location of the maxima. 
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FIG. 7: Fourier transform of the LSW for = 0.7, Vc/vp « 1.43, Vc/vs « 3, T/krVp = 2.6 x 10"*, 
same graph with charge (soUd Unes) and spin (dashed lines) wave dispersions superimposed (b). 



0.97 (a) and the 



3. Thermal effects 

On general grounds one expects that thermal effects become visible only for energies lo < T and distances x > 
At = Vs/T. Choosing e.g. T < lOK and vp « lO^m/s (a typical value for a quasi-lD organic metal for which LL 
theory should be applicablei) this implies that the spin-charge peak structure as seen in Fig. |31 will remain intact for 
distances not too far from the edge {2x//3vs ^1). We should caution the reader that the energy range for which the 
LL theory is applicable to a specific experiment may sometimes be smaller than the range depicted in Fig. [S] (cf. our 
discussion in Sec. 1). For 2x/f3vc > 1 the spin and charge waves loose their coherence before reaching the edge and 
therefore spin-charge separation is destroyed (see Fig.|^a)). 
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FIG. 8: [color online] A{uj,x;l3) for Kc = 0.7, Vc/vf ~ 1.43, Vc/vs « 3, a; = 10a, Ue = 0.97 at different temperatures. 

A more interesting issue is the fate of the asymptotic scaling behavior of the LSW near the Fermi level as the 
temperature increases. In earlier work it was found that the "uniform part" of the LSW (corresponding to the first 
term in (|34|) crosses over to a;^ scaling in both boundary (2ujx/vs <C 1) and bulk regimes {2ujx/vc 3> 1) when ujP < 1— . 
The effect was found to originate in the exponential damping of the density correlations for unequal times due to 
thermal fluctuations. By performing an expansion of the full LSW in Eq. H34|) with the small parameter luP we find 
that the power law is now modified to 



A{uj,x;f3) A + Blo + Cu^, w/3 < 1 



(36) 
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where A, B and C depend on the temperature and the distance from the boundary. We depict the crossover from 




FIG. 9: [color online] A{u;,x;(3) for T/kpVp = 2.6 x 10"^, = 0.7, v^/vp ^ 1.43, v^/vs ^ 3, = 0.97 (a) and cut at x = 200a 
(b). 

boundary scahng ^ cj'^^ ^-i)/2 for uj/3 ^ 1 to tlicrmal scaling ^ A + Buj + Cuj'^ for ujjS <C 1 in Fig.|Sl Note that one is 
able to observe this crossover for 2x/ jSvg <C 1. For 2x/ (3vc ^ 1 the boundary scahng is completely washed out. Since 
in this regime GLL{t, x, —x; P) in Eq. (|Sn|l is exponentially suppressed compared to GLL{t,x,x; P) in (j^ . it follows 
that A{uj, x; (3) ^ A + Cuj"^, and the hne shape becomes symmetric (see Fig. Eb). 

In this context it is important to point out that our finite-temperature results are strictly valid only for the case 
where an edge is modeled by an open boundary. An impurity, on the other hand, is faithfully represented by an open 
boundary condition only in the zero-temperature limit^. For finite temperatures a new energy scale Tq appears, 
characterizing the crossover from weak to strong coupling ("open boundary fixed point"), and depending on the 

impurity strength Vq. A simple RG estimate shows that To scales with Vq as To ~ V^^'"^ ^''\ with an overall scale 
factor that depends on the details of the regularization procedure^. For finite T with T ^ Tg, an open boundary 
representation of the impurity is still expected to capture the essential physics. An interesting, albeit technically 
challenging project would be to redo the calculations in this section for an impurity in the presence of the operators 
that appear away from the T — Q open boundary fixed point, allowing for a complete picture of impurity thermal 
effects. 



B. Properties of the Local Tunneling Conductance 

The local differential tunneling conductance, defined in Eq. (Q, exhibits the very same features as the LSW. The 
only difference is that the fine structure of the differential conductance is thermally smeared via the temperature 
dependence of the Fermi-Dirac distribution. This can clearly be seen by examining the energy dependence at some 
fixed distance from the boundary, as done in Fig. llUI and then comparing to the corresponding graphs for the LSW in 
Fig. As the smearing occurs on the scale of T, spin-charge separation is wiped out for temperatures T > As. The 
distance x from the boundary determines the wavelength of the oscillations in energy space, so in most cases it should 
be possible to find a range for x which shows many waves in the energy interval where LL theory applies that are 
not washed out by temperature, i.e. 1/A ^ x/wc ^ ^/T, where A is the bandwidth. Note that the inequality above 
coincides with the criterion for observing spin-charge separation in the LSW, discussed in the previous subsection. 

In wires of finite length / the energy levels become discrete and are given by a characteristic direct product of two 
spectra with uniform spacing irvc/l and ttVs/ P^^'^ ■ The standing waves of the individual levels show the corresponding 
intereference of spin and charge excitationsi^. If the temperature becomes comparable to the level spacing, a continous 
pattern of Friedel oscillations as shown above emerges. However, in order to see the predicted behavior at all, the 
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wire must always be long enough so that many levels lie within the energy interval in which the LL theory applies, 
i.e. TTVc/l ^ A. For metallic materials this translates into lengths of typically lOOnm or more. 




FIG. 10: dI/dV{y,x-l3) for T/kpVF = 2.6 x 10"^ = 0.97, with x = 50a, = 0.7, Vc/vp « 1.43, Vc/vs = \c/\s ~ 3 (a) 
and X = 100a, = 0.9, Vc/vp ^ 1.11, Vc/vs = \c/\, ~ 1.26 (b). 



When applying the results to SWCNTs similar features can be expected as has already been partially confirmed 
experimentally2Si. The dominant effect is the complicated interference pattern of the Bloch waves, which strongly 
depends on the geometry of the boundary condition and the chirality of the tubes^. Nonetheless, an enhanced 
velocity, a suppressed spectral weight and a characteristic power law of decaying oscillations are clear signatures of 
interaction effects which all have been seen in experiments^. A more complete theoretical analysis would also have 
to include the effects of backscattering, band structure, longer-range interactions and the mixing of the channels near 
the boundary which we defer to a future publication. 

V. SUMMARY 

We have derived the full finite-temperature LSW for a Luttinger liquid with an edge or impurity (modeled as an 
open boundary condition), relevant to high-precision STM measurements. We have also generalized our approach to 
the "two-channel" case that describes SWCNTs in the Luttinger liquid regime, which is qualitatively similar to the 
single- channel case. 

The LSW (determining the local differential tunneling conductance in STM measurements) exhibits a very rich 
structure as a function of temperature, distance from the impurity, and the strength of the electron interaction. 
Depending on the choice of parameters one is able to see asymmetric Fano-like line shapes and spin and charge 
oscillations. The Fano-like asymmetries are caused by an interplay of boundary and interaction effects, and, as we 
have shown, are closely linked to the Friedel oscillations in real space. Spin and charge oscillations appear due to the 
interference of propagating spin and charge waves reflected from the boundary. 

We have discussed how to consistently determine the key parameters of a Luttinger liquid (the interaction parameter 
Kc, and the spin and charge velocities) from experimental measurements of the tunneling conductance. We have also 
extensively discussed various thermal effects with focus on their influence on the Fano-like asymmetries and spin-charge 
oscillations in the LSW as a function of energy. The thermal suppression of the coherence of spin and charge waves 
makes it hard to detect interaction effects in the LSW and even more so in the local differential tunneling conductance 
(where the finite-temperature LSW gets weighted by the Fermi-Dirac distribution): For temperatures T > Vs/x (where 
Vs is the speed of the spin excitations, and x is the distance to the edge or the impurity) the characteristic power law 
behavior near the Fermi level is completely washed out and replaced by an interaction-independent analytic scaling. 
For these temperatures spin-charge separation also becomes virtually impossible to detect. 

In conclusion, our results provide guidelines for identifying and interpreting signals of electron correlations in STM 
data on SWCNTs and other one- dimensional systems, and as such should be useful in the search for realizations of 
Luttinger liquid physics. 
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